APL  PHOTONICS  1,  071301  (2016) 


Guided  acoustic  and  optical  waves  in  silicon-on-insulator 
for  Brillouin  scattering  and  optomechanics 

Christopher  J.  Sarabalis,  Jeff  T.  Hill,  and  Amir  H.  Safavi-Naeini 

Ginzton  Laboratory,  Department  of  Applied  Physics,  Stanford  University,  Stanford, 

California  94305,  USA 

(Received  29  April  2016;  accepted  10  June  2016;  published  online  1  August  2016) 


We  numerically  study  silicon  waveguides  on  silica  showing  that  it  is  possible  to 
simultaneously  guide  optical  and  acoustic  waves  in  the  technologically  important 
silicon  on  insulator  (SOI)  material  system.  Thin  waveguides,  or  fins,  exhibit  geomet¬ 
rically  softened  mechanical  modes  at  gigahertz  frequencies  with  phase  velocities 
below  the  Rayleigh  velocity  in  glass,  eliminating  acoustic  radiation  losses.  We  pro¬ 
pose  slot  waveguides  on  glass  with  telecom  optical  frequencies  and  strong  radiation 
pressure  forces  resulting  in  Brillouin  gains  on  the  order  of  500  and  50000  W  1  m  1 
for  backward  and  forward  Brillouin  scattering,  respectively.  ©  2016  Author(s).  All 
article  content,  except  where  otherwise  noted,  is  licensed  under  a  Creative 
Commons  Attribution  (CC  BY)  license  (http://creativecommons.Org/licenses/by/4.0/). 
[http://dx.doi.Org/10.1063/l.4955002] 


Stimulated  Brillouin  Scattering  (SBS)  arises  from  the  interaction  of  propagating  acoustic  and 
optical  fields.  In  many  materials  including  silicon,  Brillouin  scattering  is  the  strongest  optical 
nonlinearity.1,2  These  interactions  are  thereby  a  promising  mechanism  on  which  to  base  on-chip 
nonlinear  optical  devices  compatible  with  a  rapidly  growing  silicon  photonics  toolbox. 2-9 

While  silicon  on  insulator  (SOI)  is  a  natural  platform  for  photonics,  it  is  not  as  well-suited  for 
on-chip  phononics  and  acoustic  waveguiding.  The  high  index  contrast  between  silicon  and  silica 
glass  readily  allows  for  confinement  of  optical  fields  to  the  silicon  device  layer.  But  silicon’s  relative 
stiffness  —  and  therefore  high  sound  velocity  —  makes  guiding  acoustic  waves  difficult,  motivating 
the  use  of  soft  chalcogenide  glasses  and  partial  or  complete  releases  (removal  of  the  underlying 
glass)  to  co-localize  optics  and  mechanics.  Chalcogenide  glasses  have  exhibited  the  largest  Bril¬ 
louin  gains,  but  use  of  non-standard  materials  makes  integration  with  silicon  foundry  processes 
challenging. 111  Released  and  partially  released  silicon  structures,  while  using  standard  materials, 
require  highly  sensitive  post-fabrication  processing  steps.  Furthermore,  while  release  steps  mini¬ 
mize  mechanical  coupling  to  the  environment,  they  also  isolate  devices  from  heat-dissipating 
thermal  linkages  and  on-chip  electronics. 

The  intuition  of  “index-guided”  confinement  of  mechanical  modes  fails  for  sub-wavelength  struc¬ 
tures.  Just  as  a  beam  can  be  made  compliant  if  made  thin,  thin  silicon  strips  or  fins  on  glass  can  be 
made  to  confine  mechanics  despite  bulk  silica’s  low  sound  velocity.  This  is  related  to  how  Rayleigh 
waves  are  confined  to  a  surface  —  free  space  boundaries  affect  optics  and  mechanics  in  fundamentally 
different  ways.  For  optics,  light  propagates  faster  in  vacuum  than  in  a  homogeneous  dielectric.  In  this 
sense,  the  vacuum  is  “stiffer”  than  the  dielectric.  The  opposite  is  the  case  for  acoustic  waves;  free 
space  boundaries  soften  the  structural  response,  decreasing  the  phase  velocity  of  guided  waves. 

As  we  will  show  here,  by  creating  a  wavelength-scale  guided-wave  structure  with  carefully 
designed  free  boundaries,  we  can  in  principle  confine  both  light  and  sound  in  unreleased  silicon  on 
insulator  waveguides.  This  confinement  leads  to  low  loss  acoustic  waves  with  large  optomechan¬ 
ical  coupling.  We  present  analysis  and  simulations  showing  the  existence  of  low-loss  mechanical 
modes  confined  to  silicon  waveguides  on  a  silica  slab.  We  also  demonstrate  designs  for  struc¬ 
tures  with  simultaneously  bound,  strongly  interacting  photonic  and  phononic  modes.  These  modes 
exhibit  broadband  backward  Brillouin  scattering  without  acoustic  radiation  into  the  substrate  as 
well  as  low-loss,  forward  Brillouin  scattering.  These  results  open  the  way  for  a  new  class  of  fully 
CMOS-compatible  silicon  photonic  devices  that  utilize  mechanical  degrees  of  freedom. 


©Author(s)  2016. 
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FIG.  1.  A  SOI  slot  waveguide  with  height  h  =  340  nm,  waveguide  width  w  =  1 00  nm,  and  gap  g  =  40  nm  is  optically 
single-moded  in  the  telecom  C-band  and  supports  gigahertz  frequency  mechanical  modes  for  on-chip  backward  and  forward 
Brillouin  scattering  in  a  CMOS-compatible  architecture. 


The  interaction  between  propagating  light  and  sound  in  a  translationally  invariant  waveguide, 
as  shown  in  Figure  1,  is  well  described  by  the  coupled  mode  theory  of  SBS.  The  propagating 
electromagnetic  and  acoustic  waves  are  described  by  vector  fields 

E(x,  y,z,t)  =  ^ai(z,t)ei(x,y)elfiiZ~ia>it  (1) 

i 

and 

V(x,y,z,t)  =  b(z,t)u(x,y)e(-iK-2>--int,  (2) 

respectively.  These  waves  propagate  along  z  as  shown  in  Figure  1.  The  envelope  functions  cijiz.J) 
and  b(z,t)  capture  slow  variations  of  the  amplitude  of  the  modes  such  that  |5-a,(z,f)l  ^  \Pi<ii{z,t)\ 
and  |dfa,(z,f)|  «  |w,a;(z,f)|. 

For  there  to  be  significant  interactions  both  energy  and  momentum  must  be  conserved.  Consider 
an  incoming  optical  mode  of  frequency  a>p  and  a  mechanical  mode  of  frequency  Q.  These  modes  can 
scatter  into  a  downshifted  optical  Stokes  wave  with  frequency  a>s  -  a>p  -  Q  or  an  upshifted  anti-Stokes 
wave  with  frequency  wp  +  O.  Momentum  conservation,  or  spatial  phase-matching,  for  the  Stokes 
process  demands  /3p  =  /js  +  K.  The  process  is  commonly  referred  to  as  backward  and  forward  SBS11 
(BSBS  and  FSBS)  in  the  case  of  counter-  and  co-propagating  optical  modes,  respectively.  Addition¬ 
ally  sound  can  scatter  light  between  different  bands  of  a  waveguide,  referred  to  as  intermodal  SBS 
or  stimulated  inter-polarization  scattering  (SIPS),12  which  we  do  not  consider  here.  These  conserva¬ 
tion  laws  are  reflected  in  the  dynamics  of  the  envelope  functions  as  described  by  the  coupled  mode 
equations13 


(vpdz  +  dt)  av  =  -ig0*asb,  (3) 

(vsdz  +  d,)  as  =  -ig0apb\  (4) 

(vhdz  +  d,  +  vb(a/2  +  id))  b  =  -ig0apa*s,  (5) 

where  up,  us,  and  vb  are  the  group  velocities  of  the  pump,  stokes,  and  mechanical  wave,  a  is  the 

attenuation  constant  for  acoustic  waves,  6  is  the  phase  mismatch  K  -  ( fip  -  /3S),  and  go  is  the  op¬ 
tomechanical  coupling.  These  relations  hold  so  long  as  the  envelope  functions  vary  slowly  relative 
to  [3  and  oj  (or  K  and  O  for  mechanics).13  Assuming  the  interaction  is  weak  enough  such  that  apa* 
is  nearly  constant  over  the  mechanical  attenuation  length  or-1,  the  mechanics  mediates  an  effectively 
local  interaction  between  the  optical  modes  and  (3)-(5)  can  be  solved  for  the  steady-state  dynamics 
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of  the  optical  fields: 


3-|ap|2  =  -2Re 
dz\as\2  =  2Re  | 


Igol2 


(  t>btip(a/2  +  id) 

Igol2 

vbus(a/2  +  id) 


The  Brillouin  gain  parameter  (2Re{..})  appearing  in  these  equations  depends  on  the  normalization 
of  the  mode  amplitudes.  For  the  rest  of  this  paper,  we  assume  power  normalization  such  that  the 
calculated  gain,  T,  has  units  W  'm  1  and  is  given  by13 


r  = 


2(jjpQ.  a/ 2 

‘Pp'Ps'Pm  (a/2)2  +  d2 


/ 


2 


dA  e*  dus  u  ep 


(6) 


with  modal  powers  iP{p,s,m}. 

Silicon’s  transverse  sound  speed  c^si  =  sjE/2(\  +  v)p  =  5390  m/s  is  much  higher  than  in  glass 
where  cbg  =  3690  m/s.  The  material  properties  E,  p,  and  v  are  Young’s  modulus,  density,  and 
Poisson’s  ratio,  respectively.  Vibrations  of  a  silicon  waveguide  with  phase  velocities  above  c, ,,  will 
readily  excite  surface  and  bulk  waves  of  the  glass  substrate.  Despite  a  relatively  high  material  stiff¬ 
ness,  high  aspect-ratio  silicon  structures,  as  in  Figure  2,  can  have  geometrically  softened  modes  with 
lower  phase  velocities  than  bulk  and  surface  acoustic  waves  in  glass.  We  focus  on  strip  waveguides, 
or  fins,  in  340  nm  thick  SOI. 

Consider  the  vibrations  of  a  100  nm  wide  fin.  To  simplify  the  analysis  we  first  replace  the 
glass  substrate  with  fixed  boundary  conditions.  Solutions  to  the  equations  of  elasticity  are  computed 
numerically  in  COMSOL14  for  each  wavevector  K  by  imposing  Floquet  boundary  conditions  on 
faces  perpendicular  to  the  axis  of  propagation.  A  representative  mechanical  mode  and  the  dispersion 
of  the  clamped  beam  is  shown  (dashed)  in  Figure  2(a).  At  K  =  0  solutions  reduce  to  the  those 

of  a  cantilever  for  which15  =  J f^)2.  Here  I  =  hw3/ 12  is  the  bending  moment  of  the 

beam,  h  and  w  are  the  beam  height  and  width,  and  A  is  the  area  of  the  profile.  We  model  silicon 
as  an  isotropic  material  with  properties  E  =  165  GPa,  p  —  2330  kg/m3,  and  v  =  0.22. 15  As  K  is 
increased  the  rectangular  waveguide  exhibits  quadratic  dispersion  which  is  sensitive  to  h  and  w. 
Waveguide  geometry  becomes  insignificant  in  the  high-V  limit  where  the  mechanical  wavelength 
is  small  compared  to  w  and  h.  In  this  limit  the  phase  velocity  limits  to  cLsi  governed  only  by  bulk 
material  properties. 

We  now  consider  fins  on  glass.  The  spectrum  for  this  structure,  as  shown  in  Figure  2(a),  is 
comprised  of  flexural  modes  of  the  beam  (the  solid  line  falling  below  the  dashed  clamped  beam 


0  0.5  1  1.5  2  2.5  3  0  0.5  1  1.5  2  2.5 

K/  2n  (jxm-1)  p/2n  (nm-1) 

FIG.  2.  (a)  Mechanical  bands  of  an  h  =  340  nm,  w  =  100  nm  fin  for  clamped  (dashed)  and  on-glass  (solid)  models  contain 
intervals  of  radiation-free  modes  falling  below  the  Rayleigh  cone  (hatched).  Supermodes  of  a  slot  waveguide  (see  Figure  1) 
are  nearly  degenerate  with  the  on-glass  fin  band,  (b)  Optical  bands  of  a  ( h  =  340,  w  =  100,g  =  40)  nm  slot  waveguide  show 
a  single  TM-like  mode  (solid)  in  the  telecom  C-band.  The  fundamental  TE-like  band  (dashed)  emerges  at  higher  frequencies 
due  to  the  high  aspect  ratios  necessary  for  the  confinement  of  mechanics.  Inset  is  a  plot  of  the  electric  field. 
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dispersion  discussed  above),  and  surface  waves  at  the  glass-air  interface  (the  hatched  region)  as 
well  as  bulk  waves  in  the  glass  (also  within  the  hatched  region).  Notably  for  w  =  100  nm,  a  broad 
interval  of  the  flexural  band  lies  below  the  dispersion  of  the  surface  acoustic  waves  described  by 
the  Rayleigh  cone  with  phase  velocity  cr.  Scattering  from  this  interval  into  the  Rayleigh  waves  is 
forbidden  by  phase-matching  considerations.  Here  the  Rayleigh  wave  velocity  cr  is  related  to  the 
material  properties  of  glass  by  the  Rayleigh-Lamb  equation. 15  Modeling  the  glass  with  material 
parameters  (E  =  70  GPa.p  =  2200  kg/m3,v  =  0.17),  we  find  cr  =  3340  m/s.  The  Rayleigh  waves 
travel  slower  than  longitudinal  (5848  m/s)  and  transverse  (3690  m/s)  waves  in  glass  which  lie  above 
the  Rayleigh  line. 

As  the  width  of  a  rectangular  waveguide  is  increased,  the  flexural  band  rises  in  frequency.  At 
w  =  250  nm,  the  band  lies  completely  above  the  Rayleigh  line  and  no  protected  region  is  present. 
Owing  to  the  scale  invariance  of  the  equations  of  elasticity  and  the  dispersion  of  acoustic  waves, 
this  characteristic  aspect  ratio  250:340  depends  only  upon  material  properties  of  the  device  and 
substrate  layer  and  is  therefore  relevant  to  unreleased  SOI  devices  of  any  thickness.  For  example, 
for  220  nm  SOI,  waveguides  of  width  smaller  than  160  nm  are  required  to  achieve  non-leaky  or 
phase-protected  propagation. 

The  reflection  symmetry  of  a  fin  waveguide  suppresses  intramodal  Brillouin  scattering  via  the 
acoustic  waves  studied  above.  Under  reflection  through  the  yz  symmetry  plane,  0  V,  all  nondegen¬ 
erate  optical  and  mechanical  modes  are  either  symmetric  (+)  or  antisymmetric  (-),  i.e.,  U  (0vr) 
=  ±0AU  (r).  Consider  the  symmetry  of  the  Brillouin  gain  of  Eq.  (6).  For  intramodal  scattering, 
since  both  optical  modes,  es  and  ep,  have  the  same  symmetry,  their  product  —  and  therefore  the  op¬ 
tical  force  they  induce  —  is  symmetric.  Only  symmetric  motion  can  be  driven  by  these  symmetric, 
intramodal  optical  forces. 

The  flexural  modes  shown  in  Figure  2(a)  are  antisymmetric.  They  can  only  couple  optical 
modes  of  opposite  symmetry  and  cannot  mediate  intramodal  scattering.  In  contrast,  the  symmetric 
mechanical  supermode  of  side-by-side  fins,  i.e.,  a  slot  waveguide,  can  mediate  intramodal  scatter¬ 
ing.  Tight  confinement  of  the  fin  flexural  modes  to  the  silicon  causes  the  symmetric  and  antisym¬ 
metric  supermodes  of  the  slot  structure  to  be  nearly  degenerate  even  at  spacings  below  100  nm, 
since  very  little  acoustic  coupling  occurs  through  the  glass.  This  results  in  a  small  MHz-scale 
splitting,  not  noticeable  in  Figure  2(a). 

The  slot  waveguide  is  a  popular  photonic  structure.  The  TE-like  modes  of  slot  waveguides  are 
commonly  used  for  sensing  as  electric  field  amplitudes  are  large  in  the  gap.1'1  A  recent  study  of 
SBS  in  chalcogenide  glasses  filling  a  slot  waveguide  has  utilized  this  electric  field  enhancement 
to  achieve  large  interactions  with  index-guided  mechanical  modes.17  Slot  waveguides  have  also 
been  proposed  for  SBS  in  suspended  silicon  structures.18  Furthermore  slotted  structures  have  been 
implemented  in  photonic  crystal  fibers19  to  achieve  large  optomechanical  coupling. 

Consider  the  optical  guided  modes  for  100  nm  wide  fins  spaced  40  nm  apart  in  340  nm  thick 
silicon  shown  in  Figure  2(b).  For  these  dimensions,  the  waveguide  has  a  single,  TM-like  mode 
in  the  telecom  C-band.  The  large  splitting  of  the  fundamental  TE-like  and  TM-like  modes  of  this 
waveguide  is  a  consequence  of  the  high  aspect  ratio  necessary  to  confine  mechanics.  Additionally, 
we  expect  that  by  using  TM  guided  modes,  lower  scattering  due  to  smaller  fields  at  etched  surfaces 
will  lead  to  lower  optical  losses. 

Having  defined  the  acoustic  and  optical  modes  of  interest,  we  now  calculate  the  achievable 
Brillouin  gains  using  Eq.  (6).  The  displacement  field  causes  a  perturbation  in  the  permittivity  dus 
with  a  bulk  contribution  from  the  photoelastic  effect  duspe  and  a  boundary  contribution  from  radi¬ 
ation  pressure  <9uSrp.  For  isotropic  materials  like  silicon  and  glass,  <9t,£peu  =  -e2pS.  The  rank-4 
photoelastic  tensor  p  relates  the  strain  Sy  =  /  ((),«,  +  djU,)  to  a  change  in  e.20  This  term  is  inte¬ 
grated  over  the  bulk  of  each  dielectric.  A  second  contribution  due  to  motion  of  the  boundaries 
between  dielectrics  takes  the  form  of  an  integral  over  these  surfaces  with  integrand,21 

e*  dus rpU  ep  =  (A s  |ej_|2  -  Ae~'|d|||2)  n  •  u,  (7) 

with  the  electric  field  perpendicular  |ej2  =  e*±epx  and  electric  displacement  field  parallel  dy  to  the 
normal  of  each  boundary.  The  permittivity  differences  As  and  As  1  are  taken  from  outside  to  inside 
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in  the  direction  of  the  normal  n.  Although  in  general  both  effects  are  significant  in  wavelength-scale 
waveguides,22  radiation  pressure  at  the  silicon-air  boundary  is  an  order  of  magnitude  greater  than 
all  other  bulk  and  boundary  contributions  for  the  slot  waveguides  considered  and  the  values  of  F 
calculated  here. 

As  the  detuning  5  -  K  -  2/3  in  Eq.  (6)  is  varied,  T  sweeps  out  a  Lorentzian  with  FWHM  a 
and  height  equal  to  the  peak ,  or  phase-matched ,  Brillouin  gain  rpm  =  r|  5=().  For  each  optical  wave 
vector  fi  (Fig.  3(a))  and  corresponding  resonant  mechanical  wave  K  =  -2/3  (Fig.  3(b)),  we  compute 


.400  - 


Pl2n  ([im"1) 


FIG.  3.  For  BSBS,  resonant  optical  and  mechanical  modes  in  (a)  and  (b)  are  aligned  with  the  corresponding  peak,  or 
phase-matched,  Brillouin  gain  rpm  (c).  The  radiation  pressure  (boundary)  contribution  to  the  optomechanical  coupling  is 
an  order  of  magnitude  larger  than  the  photoelastic  effect  (body),  plotted  inset  to  (c)  on  separate  scales  for  clarity,  (d)  The  peak 
FSBS  gain  is  shown  alongside  the  872  MHz,  K  =  0  mechanical  mode.  FSBS  peak  gains  are  ~200x  greater  than  BSBS  gains, 
largely  a  consequence  of  lower  mechanical  frequency.  We  assume  a  mechanical  quality  factor  of  Q  =  1000  in  all  calculations. 
All  edges  have  a  30  nm  radius  or  fillet  for  FSBS  and  BSBS  simulations  (geometry  in  (c)).  Artifacts  of  the  simulation  have 
been  removed  from  the  curve  in  (c)  (see  supplementary  material).^ 
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the  phase-matched  Brillouin  gain  for  BSBS  (Fig.  3(c)).  This  gain  F pm  ( /3)  increases  steadily  from 
zero  at  ySmin/2;r  »  0.9  gm_1  before  turning  over  approximately  where  the  flexural  mode  enters  the 
Rayleigh  cone.  Note  that  the  gain  continues  to  increase  after  the  fin  mode  enters  the  leaky  region 
at  PI2n  >  1.4  /im  1 .  This  is  because  the  mechanical  mode  remains  localized  to  the  fin  at  larger  K 
(see  supplementary  material  for  mode  profiles).2 3 4 5 6 7 8 9  '  On  the  other  hand,  as  fi  approaches  fim]n  and  the 
optical  band  approaches  the  glass  line,  the  gain  becomes  negligible  as  the  optical  mode  spreads  into 
the  silica.  In  addition,  we  sweep  the  geometry  (see  supplementary  material23)  and  show  that  rpm 
increases  with  decreasing  w  and  g. 

As  discussed  above,  the  antisymmetry  of  the  flexural  modes  of  a  rectangular  waveguide  moti¬ 
vates  the  use  of  symmetric  supermodes  of  the  slot  waveguide  for  intramodal  SBS.  In  Figure  3(c), 
a  plot  of  the  photoelastic  and  radiation  pressure  distributions  shows  the  drawbacks  of  the  inherent 
asymmetry  of  the  flexural  modes.  These  distributions  are  roughly  antisymmetric  about  the  yz 
symmetry  plane  of  each  fin,  cancelling  in  an  optical  field  of  constant  intensity.  Despite  these  cancel¬ 
lations,  we  obtain  peak  BSBS  gains  approaching  600  W  1  m  1  that  are  comparable  to  silicon-based 
demonstrations. 

For  BSBS  a  large  interval  of  the  mechanical  band  is  forbidden  from  radiating  into  the  glass  by 
phase-matching  considerations.  In  contrast,  the  K  =  0  mechanical  mode  radiates  into  surface  acous¬ 
tic  waves  and  bulk  acoustic  modes  of  the  glass.  Due  to  impedance  mismatch,  these  loss  channels 
do  not  necessarily  prohibit  the  use  of  slot  waveguides  for  FSBS.  Simulations  of  a  w  =  100  nm, 
g  =  40  nm  slot  waveguide,  detailed  in  the  supplementary  material,23  show  that  acoustic  radiation 
limits  the  quality  factor  to  Qrd d  =  1410.  The  872  MHz  mechanical  resonance  at  K  =  0  radiates  into 
glass  modes  at  a  rate  of  620  kHz.  Current  demonstrations  of  released,  on-chip  systems  exhibit 
mechanical  quality  factors  below  this  radiative  limit.6,9 

Since  the  modal  power  is  not  well  defined  around  K  -  0,  we  express  Equation  (6)  for  FSBS 
on  resonance  as  rpm  =  y,2^.  |/  dA  e*  due  u  ep|~  and  use  this  to  compute  the  rpm(/I)  curve  in 
Figure  3(d).  The  resulting  peak  Brillouin  gain  is  two  orders  of  magnitude  larger  than  correspond¬ 
ing  values  for  BSBS  in  agreement  with  the  D.  2  scaling  of  F  from  the  mechanical  mode  energy 
normalization. 

We  have  shown  that  it  is  possible  to  make  waveguides  that  simultaneously  guide  interacting 
acoustic  and  optical  waves  in  the  technologically  important  silicon  on  insulator  platform  and  pre¬ 
sented  a  design  for  optically  single-moded,  Brillouin  active  waveguides  in  340  nm  SOI.  These 
results  open  a  significant  new  space  of  design  freedom  —  the  incorporation  of  nanomechanics 
and  acoustics  in  CMOS-compatible  photonic  architectures,  which  can  have  applications  in  optical, 
RF-microwave,  as  well  as  novel  sensing  circuits. 
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